Raman Coupling of Zeeman Sublevels in an Alkali Bose Condensate 
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We investigate amplitude and phase control of the components of the spinor order parameter of a 
87 Rb Bose-Einstein condensate. By modeling the interaction of the multilevel atomic system with a 
pair of Raman-detuned laser pulses, we show that it is possible to construct a pulse-sequence protocol 
for producing a desired state change within a single Zeeman manifold. We present several successful 
elementary tests of this protocol in both the F = 1 and F — 2 Zeeman manifolds of 87 Rb using the 
Di transitions. We describe specific features of the interaction which are important for multimode, 
spatially-varying field configurations, including the role of state-dependent, light-induced potentials. 

PACS numbers: 32.60.-K, 03.75.Mn, 42.50.Ex 



I. INTRODUCTION 

Dilute gases of alkali atoms provide a unique setting 
in which to study the properties of quantum fluids. The 
tremendous precision with which they can be controlled 
and probed makes them an ideal model system for un- 
derstanding physical phenomena important to other less 
easily studied systems with similar physical properties. 
Because alkali atoms possess internal spin degrees of free- 
dom, alkali BECs are described by a multi-component 
order parameter. This vectorial nature of the order pa- 
rameter allows for much more complex structure and dy- 
namics than occur in a single component BEC, with im- 
portant ramifications for the system's physical proper- 
ties. The work contained in this Article was motivated 
by an interest in creating and studying topological exci- 



tations of the spinor BEC such as coreless vortices 
Spin textures such as coreless vortices 0, 13, IE 
and skyrmions [E Q are important features of quan tum 
many-body systems such as liquid helium [ToL llll ] , su- 
perconductors [IH , and neutron stars [l3| , and they also 
possess important symmetry properties deeply connected 
to fundamental principles of physics 
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There have been several previous approaches to cre- 
ating nontrivial spin textures in alkali BECs. Using a 
method referred to as 'phase imprinting', the first vortex 
observed in an alkali BEC was created by coupling hy- 
pcrfine ground states using a combination of laser, radio- 
frequency, and microwave fields [IE HE [TtJ • While this 
was quite effective, the interaction time required to gen- 
erate the desired state was on the order of tens of mil- 
liseconds. In another approach, Raman transitions using 
EIT-based stopped light techniques [Hj] were used to gen- 
erate solitons and vortices. Both this approach and the 
phase imprinting experiments involve population transfer 
between different hyperfine (ground) spin states, which 
are typically separated by several GHz. This large energy 
splitting essentially freezes out spin dynamics. In a ele- 
gant and very different series of experiments, adiabatic 
ramping of an inhomogeneous magnetic field was shown 
to create spin textures within a single Zeeman manifold 



[13 . l2fj| . l2ll |22j. However, the spatial complexity of the 
spin textures that can be created using purely magnetic 
techniques is restricted by practical constraints on cre- 
ating complex field geometries and by conservation of 
angular momentum within the essentially closed system. 

In this work we explore the interaction of Raman- 
detuned laser fields with an alkali BEC, with the purpose 
of developing general all-optical techniques for controlling 
the multi-component order parameter [231 ] . We choose an 
all optical approach because state manipulation can be 
readily realized on a microsecond time scale, and we focus 
on spin states within a single hyperfine manifold because 
this provides full access to spin mixing dynamics [24l |. 
Because of the true multilevel nature of the alkali atoms, 
the atom-field interaction exhibits a remarkable degree of 
complexity which, as we will show, must be taken into ac- 
count. An interest in multilevel Raman transitions of this 
type is not new. In earlier work, for example, the cou- 
pling between three separate electronic states with Zee- 
man substructure has been studied HE HE (2?| ■ Further, 
a number of elegant experiments showing coherent con- 
trol and creation of superpositions of Zeeman sublevels 
has been conducted with atomic beams of metastable 
neon [IE [2t| [3(j. We also note that the idea of using 
a stimulated Raman process to couple multiple ground- 
state Zeeman sublevels was previously explored theoret- 
ically [3l[ and that adiabatic passage between ground 
state sublevels of a cesium atomic beam has been demon- 
strated experimentally [32| . However, to our knowledge, 
our results here represent the first detailed treatment of 
Raman coupling between sublevels of the same ground- 
state Zeeman manifold in a way that takes into account 
both the excited state and the ground state level struc- 
ture. 

One of the more interesting and useful features of the 
Raman coupling interaction investigated here is that the 
applied laser fields induce state-dependent light shifts 
which depend in a non-trivial way on the laser beam 
parameters and the atomic level structure [H, HE |35| . 
This inherent richness in the interaction provides a con- 
venient degree of flexibility in controlling the condensate, 
allowing the creation of a variety of states not accessible 
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in previous experiments. In Section |TT] we outline the 
construction of a model of the system, subject to cer- 
tain constraints which are of particular importance in 
modeling spatially inhomogeneous interactions [23[. Al- 
though the results presented here can be easily extended 
to other alkali atoms, we will focus on implementation 
in the F = 1 and F = 2 ground state sublevels of 87 Rb, 
which has a number of features that recommend it as a 
good choice for investigating BEC spinor physics. 

In Section HTT1 we show the experimental results of sev- 
eral tests of the validity of the model in 87 Rb as put 
forward in Section [TTJ These tests include demonstra- 
tions of amplitude and phase control of the components 
of the order parameter via pair- wise coupling, as well as a 
demonstration of the possibilities for simultaneous three- 
state coupling in the F = 2 manifold. The application 
of these results to the broader problem of preparing spe- 
cific spin textures with spatially-varying amplitudes and 
phases will be briefly discussed in the conclusion. 



II. MODELING THE RAMAN-COUPLED 
SPINOR BEC 

The 2 Si/ 2 ground electronic state in alkali atoms is 
composed of two hyperfine manifolds of total spin F = 
I ± 1/2, where / is the nuclear spin. These manifolds 
consist of 2^+1 sublevels, which in the absence of a mag- 
netic field are degenerate. These states are connected to 
the 2 Pi/2 electronic states by the Di transitions, which 
in the case of Rb involves photons of A = 795 nm. The 
Lande factors of the ground and excited electronic states 
differ; in the presence of a small magnetic field the tran- 
sitions between individual sublevels are no longer degen- 
erate. In constructing an appropriate model of our ex- 
perimental configuration, we consider the presence of a 
small applied magnetic bias field of up to a few tens of 
Gauss. Because the two-photon linewidth of the Raman 
coupling can easily be smaller than the nonlinear part of 
the Zeeman shift (see Section UlI Cp . it is necessary to use 
the Breit-Rabi formula to calculate the magnitude of the 
energy shifts due to the applied magnetic fields. This has 
especially important consequences when linking three or 
more states in the same Zeeman manifold. 

In both our modeling and our experiments we work 
in the low-density limit, which has two important con- 
sequences. First, the nonlinear interaction terms due to 
the mean-field energy are small enough to be neglected in 
comparison to the Raman coupling terms. Second, in a 
sufficiently low-density atomic cloud, beam propagation 
effects can be safely neglected. We also adopt a pure 
state description of the system, which is valid as long as 
the spontaneous scattering rate is negligible. This can be 
accomplished either by detuning the optical fields suffi- 
ciently far from resonance, or by using laser pulse config- 
urations which maintain the system in an adiabatically 
evolving dark state (i.e. STIRAP) [H HE H3- Spon- 
taneous heating of the condensate is undesirable in any 



case, therefore we are required to work in a regime where 
the photon scattering rate is small regardless of the state 
description. 




FIG. 1: Experimental geometry indicating the relative orien- 
tation of the Raman beams, the magnetic quantization axis, 
and gravity. The inset is a schematic indicating the linkage 
pattern characteristic of the a + - a~ coupling configuration. 

The beam interaction geometry we have chosen is 
shown in Fig [TJ The two Raman beams are copropa- 
gating, parallel to the quantization axis defined by the 
magnetic bias field. They are cr + and o~ polarized so 
that coupling occurs between states with Amp = 2. Be- 
cause the beams are collinear, the atoms' change in ki- 
netic energy due to the Raman transition is many orders 
of magnitude smaller than the condensate momentum 
distribution. This makes it possible to drop the kinetic 
energy terms from the Hamiltonian. It is also important 
to note that after the optical interaction the populations 
in different internal states have no significant relative mo- 
mentum, and remain physically overlapping after the in- 
teraction has taken place. 

Although we have limited the model and results pre- 
sented here to the a + and o~ coupling configuration, it 
should be possible to extend these results to include 7r 
transitions. This would allow all of the Zeeman sublevels 
in a manifold to be accessed, instead of the subset corre- 
sponding to Amp = 2. By using beam modes with polar- 
ization components in the direction of propagation, this 
can be accomplished while still using a collinear beam 
geo metry, i.e. by driving anomalous Raman transitions 
.'id . Experimental work exploring this intriguing possi- 
bility is in progress. 

One important feature we wish to highlight in the 
model presented here is the consideration of complex- 
valued Rabi frequencies. For interactions involving plane 
waves it is customary to constrain the Rabi frequencies 
to be real valued by factoring out an overall phase. For 
multimode configurations where the relative phase of the 
Raman beams is allowed to vary spatially, the Rabi fre- 
quencies cannot be treated as real valued over the entire 
interaction region. The experimental results presented in 
Section [IIII involve spatially uniform interactions, and so 
do not address this additional degree of freedom in the 
system, however the implications for multimode interac- 
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tions will be discussed in the conclusions. 



In the theory and experiments described in this Article, 
the detunings of the laser fields have been considered to 
lie within a range spanning about twice the excited state 
hyperfine splitting of 816 MHz. We will show that for the 
87 Rb Di transitions there is a wide range of of intensities 
and detunings for which the coherent two-photon effects 
dominate and the photon scattering rate is negligible. In 
addition, we limit our consideration to Rabi frequencies 
much smaller than the excited state hyperfine splitting, 
also noting that counter-rotating terms in the Hamilto- 
nian are typically far enough from resonance with any 
real transitions that they can be neglected. For alkali 
atoms with smaller hyperfine splittings, such as sodium, 
resonant excitation places tighter limits on the range of 
detunings that are experimentally useful. 

Bearing in mind all of these considerations, one can 
write a Hamiltonian for the system in the interaction 
picture [37| , apply the rotating wave approximation, and 
adiabatically eliminate the excited states [38| . This pro- 
cedure gives a greatly simplified effective Hamiltonian for 
the system which is valid and useful for most interaction 
configurations of experimental interest. We will investi- 
gate the application of this approach to the F = 1 and 
F = 2 ground state hyperfine manifolds of 87 Rb in the 
sections that follow. 



A. Application to the 87 Rb F = 1 Manifold 



xii++x\ c i- 



(1) 



I + and J_ are the intensities of the a + and er _ polar- 
ized beams, and 5 is the two-photon detuning defined 
below. The sub(super)scripts refer to ground (excited) 
states indicated by the labels in the state linkage diagram 
in Fig. [21 The parameters x are coupling coefficients for 
the state-dependent light shifts. Their definitions are as 
follows. 
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The quantity du t is the reduced dipole matrix clement 
for the 87 Rb Di transitions, the Cf are the Clebsch- 
Gordan coefficients for the specific transition between 
ground state i (= 1,2) and excited state a (= a, b, c, d). 
The A are detunings defined as follows. 




FIG. 2: State linkage diagram for (a + ,a~) Raman coupling 
of the F — 1 ground state manifold of 87 Rb via the Di tran- 
sitions. The magnitude of the Zeeman shift of the states by 
the small magnetic bias field is exaggerated for clarity. 
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The E n represent the atomic bare state energies, which 
depend on the magnitude of the magnetic bias field. We 
have also defined a two-photon detuning 



6 = (E 1 -E 2 )/H+(lu + -lu-) 

The parameter 77 in the Hamiltonian is a coupling coeffi- 
cient for the effective two-photon Rabi frequency. 
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The state linkages for the system composed of the 
F = 1 Zeeman manifold coupled to the Di excited states 
are shown in Fig. [2] Applying the procedure outlined 
above results in an effective two-level system, with cou- 
pling between the two indicated Zeeman sublevels de- 
scribed by a Hamiltonian with the following form. 



For completeness we have explicitly included a com- 
plex exponential term containing the difference in phase 
between the a + and o~ beams, represented as the value 
£. As mentioned in the previous section, allowing the 
Rabi frequencies to take on complex values is necessary 
for modeling interactions where the relative phase may 
vary in space or change as a function of time. 
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1. Pseudospin Representation and Fictitious Fields 

Having determined the form of the system Hamilto- 
nian, it is instructive to rewrite it in a basis of spin-1/2 
operators. 



i? = —EH er 2 + C os(£±)- o-^ + sm^ij- cr y (4) 

Here we have defined a scalar light shift a, Zeeman light 
shift 8, and two-photon Rabi freqency Q. 
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Reconstructing the Hamiltonian in the form of Eq. [4] 
gives additional insight into the expected behavior of the 
system, allowing us to visualize it as a spin-1/2 system 
subject to a fictitious electric and magnetic field [33l.l34l|. 
The scalar term involving a results in a global phase shift 
that is trivial for uniform laser fields, and has no effect on 
the internal state evolution. It should be noted, however 
that for laser fields with spatially varying intensities, the 
energy shift it represents can be physically important, 
for example it contributes to the dipole force acting on 
the atoms. The remaining terms involving spin opera- 
tors and the parameters S, f2, and £ can be conveniently 
understood as an interaction with a fictitious magnetic 
field. 



B, 



^cos(£ ± )£ + ^sin(£ ± )i+ (6 + S)k 
(j,b/H 



(8) 



If only one or the other of the laser fields is applied, 
51 = 0, and the fictitious magnetic field is oriented par- 
allel or antiparallel to the polar axis. Assuming the 
intensity-independent, two-photon detuning 6 is zero, the 
magnitude and sign are determined by that of the Zee- 
man light shift parameter d. 
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FIG. 3: Representation of the fictitious magnetic field due 
to the application of the a + laser field (left), the a~ field 
(center), and both a + and o~ fields together (right). For 
these plots, the detuning is assumed to be midway between 
the F' — 1 and F' = 2 excited state manifolds. 



solely along the polar axis. The polar angle is determined 
by the relative magnitudes of the laser field intensities J+ 
and I-. The azimuthal angle of the fictitious magnetic 
field depends on the relative phase of the fields £±. The 
relationship between the laser fields and the orientation 
of the fictitious magnetic field in the pseudospin space is 
shown in Fig. [3] for the specific case where Ar> 1 =0. 
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FIG. 4: Detuning dependence of the scalar light shift a, Zee- 
man light shift 8, and two-photon Rabi frequency fl. The 
scalar and Zeeman light shifts are expressed in terms of their 
components which depend on intensities of the a + and a~ 
polarized laser fields. (/+,/_ = 1 mW/cm 2 , B = 1.33 Gauss) 

As can been seen from Fig.[H the dependence of the pa- 
rameters a, S, and Q on detuning from the excited states 
is nontrivial. In the figure the parameters a and S are 
shown broken down into the parts which depend on the 
<7 - and <j + polarized fields. It should be noted that these 
plots are valid as long as the detuning from resonance is 
sufficiently large compared to the Rabi frequencies for 
the involved transitions. For intensities typical of our ex- 
periments (<70 mW/cm 2 ), the plots are valid except in 
a region less than 100 MHz wide around the resonances 
shown. The complexity in these curves clearly indicates 
that the detuning from the excited state must be taken 
into consideration when choosing how to act on the sys- 
tem with the laser pulses. In actual practice, we have 
found that this feature affords a degree of control over 
the interaction Hamiltonian that can be highly useful. 



2. Determining a Pulse Protocol 



When both laser fields are turned on, Q is nonzero, and 
the orientation of the fictitious magnetic field is no longer 



The motivation for the analysis of the preceeding sec- 
tion is to determine how to control the state of the sys- 
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tern, and cause it to evolve in a desired manner. Because 
the general picture of how a spin evolves in a specified 
external field is well understood, this representation is of 
great utility in understanding how to approach the prob- 
lem. A number of advanced spin manipulation techniques 
have been developed for NMR which could potentially be 
implemented in this system [39], however we focus here 
on using simple pulse configurations. 

Adiabatic passage techniques (i.e. STIRAP) can be 
used to perform coherent transfer (23[, but are not well 
suited for certain applications. For example, the large 
pulse areas required for efficient transfer make it com- 
paratively difficult to compensate for the inhomogeneous 
light shifts which arise when using beam modes with spa- 
tially varying relative intensities. We will not address the 
use of adiabatic techniqes here, focusing instead on co- 
herent population oscillation driven with square pulses. 
Specifically, we assume the use of sequences of laser pulses 
of constant intensity J_, I + , and duration r. We then de- 
termine what laser fields to apply in order to change the 
polar or azimuthal orientation of the pseudospin vector 
representing the state of the system. 

From ^ and ((7]) we see that in order to cause complete 
Rabi oscillations between the two states at a frequency 
f2, we must turn on both fields using a fixed ratio of 
intensities that makes 5 + 5 = 0. In the pseudospin space, 
this corresponds to generating a fictitious magnetic field 
which lies in the equatorial (x-y) plane, such that the 
pseudospin vector will precess in a circular path between 
1} and 2} at the antipodes. The change in the polar 
angle 9 of the pseudospin is in that case simply 



(9) 



where t is the pulse duration. It is important to note, 
however, that because the sign and magnitude of the light 
shifts depend in a complicated way on the atomic level 
structure and the detuning, it may be difficult to satisfy 
the requirement that 5 + 5 = 0. For spatially uniform 
fields, it may be possible to cancel a non-zero Zeeman 
light shift 5 over the entire area of interaction by a suit- 
able choice of the two-photon detuning i5. For a mul- 
timode interaction where the relative intensities of the 
laser fields is allowed to vary spatially, this is in general 
not experimentally feasible. As a result, in the analysis 
that follows we will assume that 5 is zero, and when cou- 
pling atoms between the states we will require the ratio 
of laser intensities to be such that 5 is zero, i.e. 
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Solving © and dTUJ) for 1+ and 7_ yields 
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FIG. 5: Detuning dependence of the intensities required to 
cause a 180° 6 rotation of the pseudospin vector describing 
the Raman-coupled \F = l,m_F = —1) and \F — l,m,F = 1) 
states, subject to the condition that 5 = 0. Negative intensi- 
ties correspond to a sign change in the direction of rotation 
of 9. Pulse duration r = 20 /is, magnetic bias field B = 1.33 
Gauss 

The relations defined in Eqns. (fTTjl and (fl"2f are plot- 
ted in Fig. |5| as a function of detuning, for typical ex- 
perimental values of pulse duration and magnetic bias 
field. It should be noted that for detunings between the 
excited states, the sense of rotation of the pseudospin 
vector reverses, which is indicated in the plot by nega- 
tive intensities. The small feature near -620 MHz occurs 
due to the fact that the Zeeman light shifts for the a + 
and a~ fields vanish there at detunings which are not 
quite identical for B ^ 0, making it impossible to satisfy 
the requirement that 5 = in the neighborhood of those 
points. 

From Fig. [5l it is apparent that using nearly equal in- 
tensities in the a + and a" fields satisfies Eq. (fTU|) over 
essentially the entire range of detunings. This is a con- 
sequence of the high degree of symmetry in the dipole 
matrix elements and state linkages. We will show later 
that in systems with less symmetry, the solutions will not 
generally be so simple. 

Having determined how to cause a 9 rotation of the 
pseudospin vector, we now consider the problem of con- 
trolling <j>, or in other words, the relative phase of the two 
states. Applying only the a + or the o~ laser fields causes 
no change in the relative population of the states, but will 
cause a change in their relative phases if the Zeeman light 
shift is nonzero. As discussed above, this effect can be 
understood in the pseudospin representation as being due 
to interaction with a fictitious magnetic field along the 
polar axis which causes the pseudospin vector to precess 
azimuthally at a rate proportional to the magnitude of 
the Zeeman light shift. The intensity required to cause 
a given rotation in <j> by applying the a + or o~ polar- 
ized beams individually for a pulse of duration t can be 
written as: 



T {e) = Vx1-X b 2 C 
+ v1Wx d 2-X b i CT 
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As with the expressions for rotation in 9, the sign of 
the rotation is determined by the sign of the light shift 
coefficients, which can be either positive or negative. For 
the F = 1 Hamiltonian under consideration in this Sec- 
tion, the light shifts due to the a + and a~ laser fields 
are opposite in sign and of nearly equal magnitude (see 
Fig. [6|). As noted above, the Zeeman light shift from 
both laser fields vanishes near -620 MHz, which makes it 
impossible to cause any <f> rotation with either laser field 
near that detuning. This is reflected in the plot by the 
obvious presence of the large singularity. It should be 
noted that only the Zeeman light shift vanishes at this 
detuning, the scalar light shift is still nonzero. 
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FIG. 6: Detuning dependence of the intensities required to 
cause a 180° 4> rotation of the pseudospin vector describing 
the Raman-coupled \F = 1, tuf = —1) and \F = l,m,F = 1} 
states. Negative intensities correspond to a sign change in 
the direction of rotation of 0. t — 20/^s, B = 1.33 Gauss 



The expressions (|1HI14I) plotted in Fig. [5] and [5] give 
a clear picture of the experimental parameter space in 
which we must work to control the system with the Ra- 
man laser pulses. There is a wide range of valid detunings 
over which we have a protocol for determining what opti- 
cal pulses will cause a desired change of the state vector. 
We will show how we have been able to apply this infor- 
mation experimentally to control the spinor wavefunction 
of an F = 1 BEC in the experimental results presented 
in section [TTTT 



B. Application to the 87 Rb F = 2 Manifold 

Having demonstrated the basic derivation of a proto- 
col for the F = 1 manifold, we turn to a more compli- 
cated case, which is to develop a useful model of Raman 
interactions within the F = 2 ground state Zeeman man- 
ifold. Several differences are immediately apparent, one 
being that there are two subgroups of the manifold which 
can be coupled by a a + , a~ beam pair: the states with 
mp = —2,0,2 and those with mp = —1,1. Because the 
rap = 2 state is the only one in which it is possible to 
produce a BEC in a magnetic trap, we will limit our at- 
tention to the mp = —2,0,2 subset, shown in Fig. [7] 
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F'= 1 




F = 2 



FIG. 7: State linkage diagram for (a + ,a~) Raman coupling 
of the F — 2 ground state manifold of 87 Rb via the Di tran- 
sitions. 



Constructing a Hamiltonian for these states in the in- 
teraction picture, making the rotating-wave approxima- 
tion, and adiabatically eliminating the excited states as 
explained above results in the following effective 3-level 
Hamiltonian. 



xtU 



77f 2 b v^O+e«± x fl_+ x fl + r,g^LT+ 







<•?= 



^23 
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As in the previous example of the F = 1 manifold, here 
we see light shift terms (xt^) m the diagonal elements 
and two-photon coupling terms (riff ) in the off-diagonal 
elements. The definition of these parameters r\ and x can 
be inferred from the sub- and superscripts by noting the 
state labels indicated in Fig. [71 and following the pattern 
of the definitions shown for the F = 1 manifold (Eqns. [21 
[3]). Several differences in form are immediately apparent, 
such as the separately defined two-photon detunings for 
the l<->2 and 2^3 transitions. 



Su = (Ei - E 2 )/h + (lo + - cj_) 
^23 = (E 2 - E 3 )/h + (uj+ - cj_) 

For B w 0, the energies of the Zeeman sublevels within 
a manifold are degenerate, and in this limit the inter- 
action terms on both sides of the"M" are equal, i.e. 
(5i2 = #i3, and 77^2 = ??23- Although it is possible to adopt 
a pseudospin- 1 representation of the system in this limit 
by expressing the Hamiltonian in terms of a suitable ba- 
sis of 3 x 3 matrices, we will consider here a different 
approach. In the presence of a magnetic field, the non- 
linear part of the Zeeman energy shifts breaks the degen- 
eracy of the two-photon transitions, i.e. #12 ^ c>23- For 
low to moderate fields (<100 G), this relative shift of the 
two transitions frequencies increases quadratically with 
the magnetic field, by 575 Hz/G 2 . For sufficiently large 
field and sufficiently small optical pulse bandwidth, the 
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l<->2 and 2<->3 transitions can be treated as separately ad- 
dressable pseudospin-1/2 subsystems with Hamiltonians 
similar to Eq. (J2}, but with the interaction parameters 
for the l<->2 transition given by 



xf) 



X?) 



i-xf 
i-xt 



Ot\2 

and the parameters for the 2 <-> 3 transition given by 



(15) 
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(16) 



As before, a is a scalar light shift, 5 is a Zeeman light 
shift, and fl is the effective two-photon Rabi frequency 
for the given transition. It is important to note that 
these pseudospin-1/2 subsystems of the F — 2 manifold 
are inherently quite asymmetric, which has a significant 
effect on the overall nature of the response of the system 
to the applied laser fields. 



F-1 



F'=2 



«23(+) 

a 23(-) 




r 20 

• 10 


h — i — r' : r-v--i 


°23(+) 

5 23(-) 

t"."-"jV. 1- 


-0.3 
=1 

''r 


-10 \ 

7 
7 

L-20 kHz ' 

F- 

r 20 

■ io j 


0.6 0.9 1.2 

A D1 (GHz) 

=2 

: 
: 
■> 

i -C.:.:.(..:.:.i.:.:.-\ 


1.2 -0.9 -0".6-. 

i 

F' 

^23 


-0.3. : 

=] 

— 1 1 


. . j 

-10 

-20 kHz ' 
F' 

r 20 
- 10/ 

— 1 — 1 — 


0.6 0.9 1.2 
Adi (GHz) 

-2 

U — i — , — i — , — i 


1 1 1 ' 1 - 

1.2 -0.9 --^^ 


-0.3 


0.3 

-10 

-20 kHz 


i 0.9 1.2 

| / A D1 (GHz) 



FIG. 8: Detuning dependence of the scalar light shift 023, 
Zeeman light shift 823, and two-photon Rabi frequency £223 
associated with the 2<->3 subspace of the F — 2 manifold. 
The scalar and Zeeman light shifts are expressed in terms of 
their components which depend on the intensities of the a + 
and (T~ polarized laser fields. (/+,/_ = 1 mW/cm 2 , B = 17 
Gauss) 

The dependence of the parameters for the 2^3 sub- 
space (Eq. I16p on single-photon detuning is shown in 



Fig. [8] The magnetic bias field assumed for these plots 
is relatively large, which results in a split in the degen- 
eracy of the two-photon transitions significantly greater 
than our typical pulse bandwidth. Although the varia- 
tion of a and f2 with detuning are generally unremarkable 
in comparison with the plots for the F = 1 system, the 
Zeeman light shifts S + and <5_ are clearly not equal and 
opposite in magnitude, as they were for the F = 1 sys- 
tem. This has important ramifications for our attempts 
to produce a protocol for controlling this system. 
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FIG. 9: Detuning dependence of the scalar light shift 012, 
Zeeman light shift 812, and two-photon Rabi frequency f2i2 
associated with the 1^2 subspace of the F = 2 manifold. 
The scalar and Zeeman light shifts are expressed in terms of 
their components which depend on the intensities of the a + 
and a~ polarized laser fields. (/+,/_ = 1 mW/cm 2 , B = 17 
Gauss) 

At zero magnetic bias field it is important to note that 
the plots of parameters (fT5)) and (fit!)) , for the l«->2 and 
2<->3 subspaces, respectively, are essentially identical if 
the identity of the a + and a~ fields is interchanged. For 
magnetic fields of more than a few Gauss, the asymme- 
try of the (non-fictitious) Zeeman shift causes significant 
differences between the two subspaces as the field is in- 
creased, particularly for detunings above resonance with 
the F' = 2 manifold. This can be observed in the differ- 
ence between Figs. HJandlHl The light shifts are much less 
dependent on magnetic field and in several ways much 
more well-behaved when the beams are detuned below 
F' = 1, which makes it a somewhat more attractive re- 
gion to work in experimentally. 

The expressions (|T5|) and (JTBJ) allow us to similarly 
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transform these subspaces into a pseudospin represen- 
tation and think in terms of fictitious fields to cause the 
system to evolve as desired. Repeating the approach used 
to generate a protocol for the F — 1 manifold produces 
a similar set of relations indicating how to cause a given 
9 or <f> rotation of the pseudospin vector within the sub- 
space. For the 2^3 transition these are: 

~-.cd I \,cd 7- 5 — ^-.cd /^.ab ~>cd 7- 

'/23 V A2 '/23 V A2 A3 



Certainly for the F = 2 manifold, choosing an appro- 
priate laser detuning is an important factor in being able 
to exert control over the state of the BEC. We have in 
general had the most reliable experimental results with 
an F = 2 BEC when operating at a detuning several 
hundred MHz below resonance with the F' = 1 manifold. 
Alternate choices of detuning are possible, but generally 
more challenging to implement experimentally. Some ex- 
amples of our approach to control of the BEC will be 
shown in the next section. 
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FIG. 10: Plot of the intensities required to cause a 180° 9 ro- 
tation (a) or 4> rotation (b) of the pseudospin vector describ- 
ing the Raman-coupled \F = 2, rtiF = 0) and IF = 2, thf = 2) 
states, subject to the condition that 823 ~ 0. Negative inten- 
sities correspond to a sign change in the direction of rotation 
of 9. Cross-hatching indicates imaginary-valued, nonphysical 
solutions. Pulse duration r = 20 /is, magnetic bias field B = 
17 Gauss 

The expressions governing the evolution of the 2<->3 
subspace in 8 are plotted as a function of detuning in 
Fig. fTUT ah The most noticeable feature of this plot is 
the large region for which the solution is not real- valued, 
and hence unphysical. Over that range of detunings S 
cannot be made zero because the sign of the Zeeman light 
shift for the cr + and cr~ laser fields is the same, and they 
cannot be made to cancel each other out. It should be 
noted that coherent transfer is still possible in this range 
of detunings, but only if the real two-photon detuning is 
adjusted such that it compensates for the Zeeman light 
shift. 

Another noticeable feature of Fig. fTUf a) is the large 
variation in the ratio of intensities required to ensure 
that 5 = 0. This can be understood as a consequence 
of the large asymmetry in the dipole matrix elements for 
the transitions involved, and the fact that varying the de- 
tuning changes the strength of coupling to the different 
excited state manifolds. 



We have performed a number of experimental tests of 
the protocols derived in the preceding sections. The ex- 
perimental procedure common to all of the experiments 
is as follows. 

A BEC of w 10 5 atoms of 87 Rb is prepared in a mag- 
netic trap, spin-polarized in either the |F = 1, nip = —1} 
or the |F = 2,mp = 2) state. This BEC is then released 
from the magnetic trap and allowed to expand for 9 ms 
to a diameter of ~ 70 fim, in the presence of a weak mag- 
netic bias field. The average density of the BEC is then 
« 5 x 10 11 atoms/cm 3 , and the mean field energy shift is 
on the order of a few tens of Hertz, which is experimen- 
tally negligible. Once this target is prepared, the Raman 
beams are pulsed on during an interaction time lasting 
5-30 /is, depending on the experiment. The diameter of 
the Raman beams is large enough (1.7 mm) to ensure 
that the beam intensity is essentially uniform over the 
area of the BEC, with AI/I avg < 0.01. The magnetic 
bias field over this region is uniform to better than 1 mG 
at the maximum applied field strength of 17 G. 

After the laser interaction, we determine the final dis- 
tribution of population among the Zeeman sublevels by 
applying a strong transverse magnetic field gradient of 
«300 G/cm for 1 ms, which causes a Stern-Gerlach sep- 
aration of the different spin components during a sub- 
sequent 20 ms time-of-flight [19] . The field gradient is 
uniform to within ±10% over the condensate, which is 
sufficient to avoid significant distortion due to magnetic 
lcnsing during the separation process. Once separated, 
the clouds are imaged using typical absorption imaging 
techniques. The population in each Zeeman sublevel is 
quantified by integrating the measured density over the 
appropriate regions of the image, correcting for variations 
in background level. The ratio of population in the target 
state to the total population can then be calculated and 
compared to theoretical predictions. The details of the 
individual experimental tests we conducted are explained 
below. 



A. Rabi Oscillations 

The primary test of the model and the protocols we 
have derived is to see whether complete population oscil- 
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lation occurs at the frequency predicted for the specified 
intensities of the tr + and cr~ fields. Although the typical 
method of observing Rabi oscillations is to apply a field 
and vary the duration of interactions, this approach is 
not ideal for testing of our model of the system. 

The approach we chose to use to observe Rabi oscilla- 
tions in our system was to hold the pulse duration con- 
stant and increase the intensities in the ratio indicated 
by the equations, observing the resulting change in the 
target state population. Any mismatch between theory 
and experiment would result in the system being driven 
increasingly out of resonance at higher pulse intensities 
by the increasing light shifts. Another motivation to per- 
form such a test is that the model was developed with 
the intention to apply it to laser modes with spatially 
varying intensity. The performance of the model in this 
respect is therefore an important metric of its utility. 




0.0 . , . , . , — . r- 

10 20 30 40 50 



/IX (mW/cm 2 ) 

FIG. 11: Population in the |1, 1) state after a 20 fis Raman 
laser pulse as a function of square root of the product of the 
beam intensities. The vertical error bars indicate the scatter 
range in the measurements, and the horizontal error bars indi- 
cate the uncertainty in the beam intensity due to laser power 
fluctuations. The grey line is the theoretically predicted in- 
tensity dependence. 

The results of an experiment demonstrating Rabi os- 
cillations between the |1, — 1) and 1, 1) states are shown 
in Fig. [11] In this experiment, we used optical pulses 
of 20 /is duration and intensities of up to 55 mW/cm 2 . 
The beams were detuned midway between the F' = 1 
and F' = 2 transitions, i.e. at Ab 1 = on the plots 
of section [H] At this detuning, the ratio of 1+ /I_ which 
makes S = is 1.01. The magnetic bias field applied to 
the system was 1.33 Gauss, resulting in a ground-state 
Zeeman splitting of 0.93 MHz. 

The plot shows the fraction of the population trans- 
ferred to the |1, 1} state as a function of the square root 
of the product of the beam intensities, which is propor- 
tional to Cl. The light grey line is the theoretical predic- 
tion from the model. Each data point is the average of 
several runs of the experiment, with the error bars indi- 
cating the full range of scatter of the individual runs. Due 



to beam intensity variations, the vertical scatter is more 
pronounced in the regions where the function is steeply 
varying, however the model and the data are in excellent 
agreement. 
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FIG. 12: Population in the |2,0) state after a 20 fis Raman 
laser pulse as a function of square root of the product of the 
beam intensities. The vertical error bars indicate the scat- 
ter range in the measurements, and the horizontal error bars 
indicate the uncertainty in the beam intensity due to laser 
power fluctuations 

The results from a similar experiment demonstrating 
coherent oscillations between the |2, 2) state and the |2, 0) 
state are shown in Fig. [T2l The magnetic bias field ap- 
plied in this experiment was 17 Gauss, resulting in a 
ground-state Zeeman splitting of 12 MHz, and a differ- 
ence between the two-photon resonances of S12 — ^23 — 
166 kHz. This splitting was large enough that with laser 
pulses of 20 /is duration (l/r=50 kHz), we were able to 
address the |2, 2) <-> |2,0) transition separately from the 
|2, -2) <-> |2,0) transition. 

For this experiment, was set to -0.8 GHz, below 
resonance with the F' = 1 excited state manifold. We 
maintained the ratio of laser intensities at I+/I- = 2.40, 
which is the value that the model predicts will result in 
S = 0. It is worthwhile to note that there is an upper limit 
to the intensities which can be used in these experiments, 
because of the necessity of avoiding incoherent exitation 
of the atoms. For the detunings and pulse durations used 
in these experiments, we typically only begin to see sig- 
nificant signs of incoherent excitation at intensities well 
above 200 mW/cm 2 . 

In both Fig. [TT] and [T2J the experimental results are 
in excellent agreement with the theoretical predictions. 
We emphasize that the same calibration values for the 
absolute intensity have been used for all the data sets 
presented in this Article, and there are no free parameters 
in any of the theoretical curves. We do note that there 
is some decrease in the amplitude of the oscillation at 
higher intensities. This could be due either to technical 
factors such as fluctuations in the beam intensities and 
drift in the magnetic bias field, or possibly small physical 
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influences neglected in our current model, such as beam 
propagation effects. If the cause is one of the latter, it 
would prove to be an interesting subject for future study. 

B. Phase Control in the F — 1 Manifold 

After having established that we can cause several 
complete Rabi oscillations between two Zeeman sublevels 
by applying an appropriately configured Raman pulse 
pair, we can approach more sophisticated problems in 
system control. A more demanding test of the protocols 
derived in Section[TT]is to attempt to demonstrate control 
of the relative phase of the states being coupled. 




= 7t/2 cp = 9 = 7t/2 



b) 




12) 2) [2) 

9 = nil 9 = it 9 = nil 



FIG. 13: Applying two successive 9 — n/2 pulses with no 
phase shift results in complete transfer to the final state (a). 
Turning on only one of the Raman lasers causes a phase shift, 
rotating the pseudospin vector in (j>. For a (j> — ir rotation (b), 
the population returns to the initial state. 

One simple way of showing phase control in these 
pscudospin-1/2 systems is to perform a Ramsey- fringe 
type experiment, using a three step pulse sequence as 
depicted in Fig. Q1J] The first step in this sequence is to 
apply a a + -a~ pulse pair which causes the pseudospin 
vector to precess 90° in the 9 direction (tt/2 effective 
pulse area). The second step is to apply only one or the 
other of the laser fields, which causes the pseudospin to 
precess in <fi by an amount which depends on the mag- 
nitude of the Zeeman light shift and the duration of the 
pulse. The last step is to again apply a 9 = tt/2 pulse. 
If neither of the fields is turned on in between the two 
9 = tt/2 pulses, after the final step all the population 
should end up in the target state. (Fig. [l~3a ) If a small 
phase shift is applied during the second step, the last op- 
eration will not move all of the population into the final 
state. If the 4> phase shift is 180° (ir), the action of the 
final pulse will be to actually return all of the population 
to the initial state. (Fig. fl"3b) 

The results of an experiment implementing the pulse 
sequence of Fig. [T3] are shown in Fig. Q31 For this ex- 
periment, we used the |1, — 1) <-> |1, 1) system, with the 
laser detunings set to Ad 1 = as in the coherent oscil- 




Phase Shift Expected for 5 = 

FIG. 14: Experimental results of applying the 9 — n/2, cj>, 
9 — 7r/2 pulse sequence described above to the |1, — 1) <-» 
jl, 1) pseudospin system. Application of a <j> changing pulse 
is shown to cause the expected oscillation in the final state 
population. The grey line is a theoretical prediction taking 
into account a nonzero S caused by lowering the a- beam 
intensity by 5%. 

lation experiment of Fig. [TT] above. The 9 — tt/2 pulses 
were 5 fj,s in duration, separated by 10 [is. After the first 
9 = tt/2 pulse, we left on either the <j + or the a~ laser 
field for a duration of 0-10 /is in order to cause the desired 
phase shift. For the data shown, the ratio of I+/I- was 
offset by 5% from the values which would make 8 = 0, in 
order to make the sign difference in the rotation caused 
by the a + and a~ beam in this configuration explicitly 
apparent. The theoretical curve in the figure has been 
shifted accordingly, and is in reasonable agreement with 
the experimental data. This confirms that phase control 
of the system is indeed possible, and occurs as predicted 
by the model. 



C. Two-photon Lineshape in the F — 2 Manifold 

In the previous Sections we have focused on Raman- 
coupling of two selected states out of a Zeeman manifold. 
The distinction made in Section lll Bl between treating the 
|2,2), 1 2, 0) , and |2, — 2) states as either a pseudospin-1 
system at zero magnetic field, or a pair of pseudospin- 
1/2 systems in the presence of a large magnetic field is, 
of course, an oversimplification. At moderate magnetic 
fields, using pulses of sufficient bandwidth, this trio of 
states in the F = 2 manifold can be coupled simulta- 
neously even though the two-photon transitions are not 
quite degenerate. This allows for the intriguing possibil- 
ity of using a single laser pulse pair to distribute pop- 
ulation into all three states, where the fraction of the 
population in each state and their relative phase can be 
controlled by varying the splitting between S12 and S23 
with the magnetic field, and adjusting the two-photon 
detuning with respect to these resonances in a manner 
determined by numerical modeling of the system. 

The results of an experiment demonstrating the two- 
photon detuning dependence of the fractional transfer to 
the |2,0) and |2, — 2) states are shown in Fig. [TSl The 



11 



o 

3 i 

Oh 



0.5 
0.4 
0.3 







• 


— 2, Theor. 








— 2,-2 Theor. 








• 2, Exper. 








□ 2,-2 Exper. 


■ ^-0 □ C^tK 


Z i • 

^? J \* I 

\0 7 \ P 

V-6 □ □ □ D y 


\ / # v 

V / • \ 



-0.2 



-0.1 



0.0 
5 (MHz) 



0.1 



0.2 



FIG. 15: Fractional population transferred to |2,0) and 
2,-2) from a BEC initially in the |2,2) state by a 20 ps 
Raman laser pulse for different choices of the two-photon de- 
tuning S. The magnitude and ratio of the intensities were set 
for an effective 6 = 1.5-7T pulse area on two-photon resonance 
with the 1 2, 2) <-» 1 2, 0) transition, with the detuning of the 
beams set at = -800 MHz. The applied magnetic bias 

field was 17 Gauss, for which <5i2 — 523 = 166 kHz. 
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applied magnetic field used for that data set is still rel- 
atively large, such that the energy separation of the two 
photon resonances is several times the pulse bandwidth. 
The intensities used were I + = 30.6 mW/cm 2 and I_= 
12.8 mW/cm 2 , which are appropriate to make 623 = 0, 
with an effective two-photon pulse area of 1.5 it for res- 
onance with that transition. It should be noted that in 
this plot the maximum vertical scale for the transfer frac- 
tion is consequently 0.5, not 1. With the system in this 
configuration, two-photon detuning was varied about res- 
onance with the |2,2) to |2,0) transition by ±200 MHz, 
and the fractional transfer of population to the |2, 0) and 

2,-2} states was recorded, taking the average of sev- 
eral experimental runs. The oscillatory behavior of the 
transfer to |2,0) appears with a periodicity close to that 
expected from the pulse bandwidth. As the two photon 
detuning is moved below -100 MHz, however the lasers 
begin to come into resonance with the transition to the 

2,-2} state, and population begins to appear there as 
well. The solid curves appearing behind the data points 
are from a numerical model of the three-state dynamics 
in this configuration, and can be seen to be in excellent 
agreement with the experimental results. 

The general dependence of this three-state coupling on 
beam intensity adds another degree of complexity to the 
picture, as can be seen in the numerical predictions plot- 
ted in Fig. [TBJ The upper plot (a) indicates the fractional 
transfer to the |2,0) state as a function of two-photon 
detuning in the broad axis, and as a function of the ef- 
fective pulse area when resonant with the 2, 2) <-> |2, 0} 
transition along the other axis. The intensity ratio and 
magnetic field are the same as those used in the results 
of Fig. rrS] The lower plot (b) shows the transfer fraction 
to the |2, —2) state for the same experimental configura- 
tion. The theoretical curves shown in Fig. [15] are a cross 
section taken from the data of Fig. [TBI as indicated by 
the white lines on the surfaces. 

Several interesting features are apparent in Fig. 1161 



FIG. 16: Numerical prediction of the population transferred 
to |2,0) (a) and |2, -2} (b) from a BEC initially in the |2, 2} 
state, as a function of two-photon detuning and effective pulse 
area. The ratio of the intensities was set to eliminate 5 for the 
|2, 2) <-> |2,0) transition. A Dl = -800 MHz, B = 17 Gauss, 
for which S12 — 523 = 166 kHz. The theoretical curves shown 
in Fig. [15] are lineouts from these plots at the indicated loca- 
tion of 6 = 1.5-7T The highest peaks corresponds to complete 
population transfer. 



the most striking of which is that for high enough inten- 
sity and a two photon detuning of about <5 12 =-160kHz, 
all the population bypasses the |2, 0} state and ends up 
in 2,— 2). This resonance ridge is noticeably narrower 
than the peaks seen in Fig. [TBJ a, and clearly moves off 
toward lower two-photon detuning as the beam intensi- 
ties are increased. These characteristics are a result of 
the increasing uncompensated Zeeman light shift of the 
2,0) <-> |2, — 2) subsystem. The ability to place popula- 
tion in all three selected states in a tunable, determinis- 
tic manner is an extremely useful feature of the Raman- 
coupling system, allowing for the creation of more com- 
plex multi-component spin textures. 



IV. CONCLUSION 

We have shown that the Di transitions in 87 Rb can 
be used to coherently control the amplitude and phase 
of selected Zeeman sublevels in both the F — 1 and 
F = 2 ground state manifolds. The model we have put 
forward accurately predicts the dynamics of the system 
and has allowed us to create simple protocols for pro- 
ducing a desired change in the atomic state. We have 
demonstrated successful implementation of these proto- 
cols both for single-pulse-pair operations, and for se- 
quences of pulses. The results presented here are specif- 
ically applied to an untrapped, freely expanding BEC, 
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however, this approach can be adapted to model Raman- 
coupling of a BEC in a trapping potential, especially if 
the pulse bandwidth is large compared to the frequency 
shifts due to the trapping potential and mean field. 

Although we have derived and tested the model pre- 
sented here in the plane- wave limit, it can be readily 
extended to describe interaction geometries employing 
laser fields with spatially- varying intensities and phases. 
For example, it is possible to create vortices in the BEC 
order parameter by using optical vortex beams with an 
azimuthal phase winding, such as the Laguerre-Gaussian 
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